Online Supplementary Material: 
Low temperature vortex liquid in La2 xSrxCu04 
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I. INTRODUCTION 

We report here the supplementary material for Ref . [l[ . 
Using torque magnetometry, we measured the magneti- 
zation anisotropy of 7 crystals of La2_xSrxCu04 (LSCO), 
which are labeled as 03 (with x = 0.030), 04 (0.040) 05 
(0.050), 055 (0.055), 06 (0.060), 07 (0.070) and 09 (0.090). 

In optimally-doped cuprates, the bulk susceptibility is 
dominated by the paramagnetic van-Vleck orbital term 
Xorb which has a significant anisotropy {x"^^ > Xah^) 
that changes weakly with T (subscripts c and ab identify 
quantities measured with H||c and H ± c, respectively). 
Moreover, in the lightly-doped regime, the paramagnetic 
spin susceptibilities Xc ^^'^ Xab become significantly large 
below the interval 40-60 K. However, the spin susceptibil- 
ity is very nearly isotropic (except below 10 K where its 
anisotropy becomes measurable). Against the large or- 
bital and spin terms, the weak diamagnetic signal is very 
difficult to resolve using standard bulk magnetometry in 
lightly-doped cuprates. By contrast, torque magnetom- 
etry selectively detects the orbital diamagnetism gener- 
ated by supercurrents confined to the Cu02 layers while 
ignoring the large spin contribution when it is isotropic. 
The orbital van-Vleck contribution is also detected, but 
as a "background" that is H-\meai to intense fields and 
only mildly T dependent. 



II. EXPERIMENTAL DETAILS 

Each crystal, of nominal size 2x1x0.35 mm'^, was 
glued to the tip of the cantilever with its c-axis at an angle 
(f> 15° to H. The torque r = m x B leads to a flexing 
of the cantilever which is detected capacitively (m is the 
sample's magnetic moment and B = /xo(H -I- M) with /io 
the vacuum permeability). We resolve m 5 x 10^^ emu 
at 10 T (I emu = 10^"^ Am^). The torque measurements 
were performed in-house in fields up to 14 T (down to 
T — 0.35 K). High-field measurements to either 33 or 
45 T were carried out at the National High Magnetic 
Field Laboratory (NHMFL), Tallahassee. Bulk measure- 
ments of all samples were performed to 2 K in a SQUID 
magnetometer (with resolution 10^^ emu) to calibrate 
the torque cantilever. SQUID magnetometry was also 
used to observe flux expulsion by measuring the Meiss- 
ner curves with H = 10 Oe after zero-field cooling. The 
curves in Fig. [T] reveal full flux expulsion in sample 07, 
but only partial expulsion at the lowest T (2 K) in 06 
and 055. No Meissner signal is observed in sample 05 



down to 2 K. The Meissner effect requires the existence 
of long-range phase stiffness and coherence. 

It is convenient to express the torque as an effective 



observed magnetization [5|, l6| Mo. 
_Bsin0 (we take z||c). We have 



t/B^V, with Bj; 



-Ax°'-''(T)iJ„ 
(1) 

where Md is the diamagnetic term of interest. Because 
the normal-state resistivity anisotropy is extremely large 
in LSCO for x <0.10 [pd Pab = 6,000-8,000 below 40 
K [4]), we may assume that the supercurrents are pre- 
dominantly in-plane. Hence the in-plane component 
Md^y is negligible especially at fields above the melting 
field of the vortex solid Hm- The term AA/^ from the 
spin response is derived in Sec. IIVI The anisotropy of 
the van Vleck orbital susceptibility /S.x°^^ is the differ- 
ence Ax°''^ — x°J^ — x'ab (hereafter, we write H for Hz). 
In LSCO, Xc ^ 0.055". 

As a check, we have compared the "raw" data Mobs/ H 
in sample 03 against the SQUID magnetometry measure- 
ments of Lavrov et al. (Fig. HD. They have reported 
detailed measurements of the total bulk susceptibilities 
Xa, Xb and Xc in a very large LSCO crystal [x — 0.03) 
with a volume ~10 times larger than in our samples. 
In Fig. [5] we plot the T dependence of the anisotropy 
= Xc — \iXa + Xb)- There is good agreement with 
our corresponding quantity Mobs/H measured at 2 and 
5 T above 40 K. 



III. MAGNETIZATION RESULTS 

In the high-field experiments at the NHMFL, all sam- 
ples (03-09) were investigated to a peak field of 33 T (us- 
ing Bitter magnets). Further measurements were done on 
Sample 06 in the hybrid magnet in the field range 11-45 
T. In all samples, the high-field data were supplemented 
by extensive "low-field" in-house results (up to 14 T), 
which have improved resolution below ~5 T. 

Figure [3] shows the general trend of the "raw" observed 
magnetization Mobs vs. H in Samples 03-07 using the 
low-field data sets. Above the onset temperature for the 
appearance of vortex-Nernst signal Tonset, (~ 30, 40, 60 
K in samples 04, 05 and 06, respectively), the _ff-linear 
term Ax°^^H dominates to produce a fan-like pattern. 
Below Tonset, the rapid growth of Md{T,H) which is 
strongly nonlinear in H produces a noticeable downward 
deviation from the fan pattern. This is strikingly evident 
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in 055, 06 and 07 but is also seen in 03, 04 and 05. When 
T decreases below T^, the samples with x > (055, 06 
and 07) display a large diamagnetism even in the limit 
H ^ 0. This requires the vortex solid to be stable. In 
samples with x < Xc (03-05), however, the diamagnetic 
signal does not grow significantly even with cooling to 
0.35 K. 

Susceptibility The 3 distinct contributions (Eq. [1]) to 
the observed magnetization Mobs are actually readily ap- 
parent in the raw magnetization data when the field is 
extended to 33 Tesla. In a linear plot of Mq^s vs. H 
(as in Fig. [3^), the terms Md and AMg tend to be 
dwarfed by the large van- Vleck term Ax°''^H, and hard 
to make out. On the other hand, by simply dividing the 
raw magnetization by H to form the observed suscepti- 
bility Xobs = Mobs/H, we can make all 3 terms apparent. 
Figure|4]shows the field dependence of Xobs {T, H) in Sam- 
ples 03-055 at selected T [we display Xobs{T,H) divided 
by its curve at 40 K]. First, we look at Panel d (Sam- 
ple 055). Cooling from 80 K to 35 K, we observe that 
Xobs is H independent with a weak T dependence. This 
high-T behavior of Xobs leads to the "fan-like" pattern 
identified with the van- Vleck magnetization Ax°^^H in 
Fig. [3l The negative contribution of Md/H, first resolved 
near 50 K, grows rapidly in magnitude as T decreases to 
5 K. Moreover, while Md/H is very large in low fields, 
it is suppressed when H exceeds 33 T, consistent with 
field suppression of the pair condensate. Sample 05 in 
Panel c shows a similar pattern of behavior except that, 
below 10 K, a new positive, strongly iJ-dependent con- 
tribution emerges to lift the curves upwards (compare 
curves at 5 and 10 K). Proceeding to Samples 04 and 
03 (Panels b and a, respectively), we see that the new 
paramagnetic term becomes steadily larger with deceas- 
ing X. This term - associated with the spin contribution 
AMs/H - is clearly distinguishable from the van- Vleck 
and diamagnetic terms. In spite of the large spin term 
in 03 and 04, the diamagnetic term Md/H remains quite 
robust up to fields of 20-30 T. 

Removal of van Vleck term As explained in the main 
text [1], it is best to subtract altogether the large "back- 
ground" van Vleck term Ax°^^ to analyze the diamag- 
netic and spin terms accurately. The broad interval of T 
in the low-field data set allows Ax°^^{T) to be measured 
accurately. As an example, Fig. [5] displays the T depen- 
dence of Mobs (measured at 5 T) in sample 06. When 
T > 50 K, Mobs is dominated by Ax°''^H which has the 
weak T-linear dependence 

Ax'"-\T) ^a(^1 + ^^ (2) 

with A ~ 1.05x10"^ and Tq - 700 K (dashed line). The 
temperature scale Tq is independent of x to our reso- 
lution, while the intercept parameter A shows a very 
weak x dependence (both are associated with the van 
Vleck background). However, Tonset decreases signifi- 
cantly, nominally linearly with falls below 0.07. 
The change between and 40 K of this background term 



is < 6%. Using Eq. [2l we remove the van Vleck term to 
obtain the curves of M'^^^iT, H), as explained in Ref. [l|. 

As anticipated in Fig. 21 the background-subtracted 
curves Af^^^ reveal an interesting competition between 
the diamagnetic term Md{T,H) and the spin term 
AMsiT,H) as x changes from 0.03 to 0.06 (Fig. ^. In 
sample 03, Md though small is still strong enough to pull 
^'^'obs negative below --10 T. Above 10 T, Md is rapidly 
suppressed to zero at 24 T, leaving the spin term which is 
nominally H independent. Turning to Panel d for sample 
06, we see the same competition between the 2 terms, ex- 
cept that Md is now much larger, and a larger depairing 
field (~48 T) is needed to suppress it to zero. Also the 
saturation value of AMs is only i as large. The behav- 
iors in 04, 05 and 055 are intermediate between these 2 
extremes. 

Paramagnetic spin-^ moments The overall behavior of 
the spin term in high fields, together with the oscillatory 
behavior in weak H (see Fig. Ic of Ref. [Ij]) suggests 
that the spin term arises from spin-^ local moments that 
are nearly non-interacting but have anisotropic (7-factors 
{9ab < <?c) that are weakly T-dependent. As derived in 
Sec. llVi the spin contribution to Mobs is 

AMsiT,H) = ViT)ta.nh[Pg^fiBB/2], (3) 

with /is the Bohr magneton and /? = l/fcsT. The pref- 
actor is {Ng is the spin population) 

nr) = ^%^cos0. (4) 

V 2^0 
with the effective g-factor 

50 = V {Qc cos (f) + {gab sin <j)) 2 . (5) 

With 50 ^ gc fixed at 2.1, the only adjustable parameter 
at each T is "P. 

Using Eq. [31 we may subtract the contribution of the 
spin term from M^j^^ to obtain the diamagnetic term Md 
at each T. The curves of both AMg and Md in sample 03 
are displayed in Fig. [7l The profile of Md resembles the 
"tilted-hill" profile observed for the vortex-Nernst signal 
observed at T > Tc in both underdoped LSCO (at large 
doping) and in Bi-based cuprates. The curves here are 
all associated with the vortex liquid (the vortex solid is 
not observed in 03 down to 0.35 K). 

We remark that the diamagnetic term Md is readily 
apparent in the raw data of Xobs vs. H in all samples. 
Moreover, it remains robust up to high fields of 25-45 T 
(Fig. m, so the evidence for the vortex- liquid diamag- 
netism in samples 03-06 does not depend on our back- 
ground subtraction or the form of Eq. [31 The merit of 
Eq. [31 is that it allows us to understand the details of the 
low-T oscillatory behavior of M^^^ shown in Fig. Ic of 
Ref. [H, and to extract Md{T,H) with good accuracy. 

In Fig. [71 the paramagnetic curves are the best fits us- 
ing Eq. [31with the sole adjustable parameter V{T) which 
is plotted in the inset. Below Tsg ~ 2.5 K in this sample. 
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the curves of M^^^, vs. H become slightly hysteretic re- 
flecting the onset of spin-glass behavior. The spin-glass 
hysteresis loops are clockwise in contrast with the anti- 
clockwise sense of the loops in the vortex solid (we will 
discuss the spin-glass observations elsewhere). 

IV. ANISOTROPIC SPIN TERM 

The Hamiltonian for an s = ^ spin with anisotropic 
g-factor in a field B (lying in the x-z plane at an angle (j) 
to z||c) is 

H = -]^liBB{gab<yxSa^4> + gcf^z cos (j)), (6) 

where Ux and are the Pauli matrices and B ~ |B|. 
Because of the anisotropy, the tilt-angle 6* of differs 
from (f). The eigen-spinors of Eq. [6] are 

w=(ri)' i-)=(;:^f) 

with 9 defined by 

cos0 = — coS(^, sin0=— sin0, (8) 
94> 94> 

and as given in Eq. [5] With Eq. [71 the matrix ele- 
ments are (+|o'z|-|-) = cos6' and (+|cra;|+} = sin 6* etc. 

For a collection of Ns spins in a volume in a bath at 
temperature T, the thermal-averaged magnetization is 

{Ms) = ^2V^ [gc cos 6>z + gab sin 6>x] iaah[(jg^^iBB /2\. 

(9) 

This contributes the torque signal 

Ts = NstiB^^^^^Bcos(j)siii(btaiih[f3g^^iBB/2l 

(10) 



from which we obtain the results in Eqs. [SHSl 



V. HYSTERETIC CURVES AND VORTEX 
AVALANCHES 

A way to see the rapid suppression of the vortex solid 
state as x — > is to compare the hysteretic curves of 
Samples 06 and 055 at low temperatures (Fig. [5]). In 
going from 06 to 055 at T = 0.35 K, the width of the 
hysteretic loops at a fixed field {H = 1 T) shrinks by 
~100. Between 07 and 06, the loops shrink by another 
factor of 10 (see Fig. [3]). The 1000- fold decrease between 
07 and 055 provides evidence for a critical transition at 
Xc- The opposing view of a distribution of superconduct- 
ing islands implies a gradual decrease of the hysteretic 
loops that is inconsistent with our observations. 

At T = 0.35 K, vortex avalanches can be triggered in 
the vortex solid {H < Hirr{T)) by sweeping the field at a 
rate of 3-5 T/min. As shown in Fig. [51 the jumps occur 
during the field sweep- up portion of the hysteretic loops. 
By decreasing the sweep rate by a factor of 10, we can 
eliminate the jumps altogether. Vortex avalanches occur 
when the field sweep is too rapid to allow the inserted 
vortices to equilibrate in the solid phase. They involve 
collective depinning of a large fraction of the vortices in 
the crystal in a very short time interval, and are generally 
considered to be direct evidence for the existence of the 
vortex solid. 

The research at Princeton was supported by the Na- 
tional Science Foundation (NSF) through a MRSEC 
grant DMR 0213706. Research at CRIEPI was supported 
by a Grant-in- Aid for Science from the Japan Society for 
the Promotion of Science. The high field measurements 
were performed in the National High Magnetic Field Lab. 
Tallahassee, which is supported by NSF, the Department 
of Energy and the State of Florida. 
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FIG. 1: The Meissner curves measured in weak field {H— 
10 Oe) in samples 05, 055, 06 and 07 measured after zero- 
field cooling. The Meissner transitions are observed with mid- 
point temperature Tc = 0, 0.5, 5, and 12 K in 05, 055, 06 and 
07, respectively. Full flux expulsion (with demagnetization 
factor Nc ^ 2.2) corresponds to Md = -1,900 A/m. It is 
instructive to compare these ZFC curves with the high-field 
curves in Fig. |31 which are all in the reversible regime (except 
below 5 K in 07). 
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FIG. 2: Comparison of the total Mobs/H measured a.t H = 
2 and 5 T in sample 03 with the susceptibility anisotropy 
Ax inferred from the bulk measurements of Lavrov et al. [31 • 
In the latter, Ax is defined as Xc — Xab where the in-plane 
susceptibility Xab = |(Xa + Xb)- The results of Lavrov et 
al.were measured by SQUID magnetometry on a crystal with 
X = 0.03 (and ~10 times larger than 03 in volume) at the fixed 
field H = 5 kOe. Above 40 K, the agreement is quite good 
(the slight offset may come from absolute calibration of the 
torque cantilever). Below 40 K, the ratio Mobs/H is strongly 
H dependent because of increased contributions from Md and 
AM,. 
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FIG. 3: The total (observed) magnetization Mobs in the 6 LSCO crystals 03, 04, 05, 055, 06 and 07 with doping x = 0.03, 0.04, 
0.05, 0.055, 0.06 and 0.07, respectively. The Meissner transition midpoint Tc ~ 0.5, 5 and 12 K in 055, 06 and 07. The range of 
T in each sample (4-100 K) are similar to the ones indicated for sample 07. All curves share the same scales for the Alots and 
H axes, but the origins have been shifted vertically for clarity. The background paramagnetic term Ax°' is strictly //-linear 
and weakly T dependent (nominally the same in all samples). Below 30 K, the growth of the diamagnetic term Md, and the 
spin anisotropy AMs, become increasing evident from 03 to 07. The horizontal bold curve plotted with sample 03 shows the 
signal of the unloaded cantilever amplified by 10. All curves are in the reversible regime, except in sample 07 where hysteresis 
appears below 5 K. See Ref. [l[ for discussion of irreversibility. 
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FIG. 4: The H dependence of the observed susceptibihty Xoba = Moba/ H at selected T in Samples 03, 04, 05 and 055 (Panels 
a-d, respectively). Xobs{T,H) is plotted divided by its value at 40 K. In each panel Xobs is H independent and only weakly 
T dependent above ^35 K consistent with the van-Vleck term. The strongly //-dependent diamagnetic term Md/H grows 
rapidly between 20 and 5 K. Finally, at low T, a third positive contribution - the paramagnetic spin term AMs/H - emerges 
to counter the diamagnetic term in intense fields. At 5 K, the spin term is barely visible in 055, but quite large in Samples 03 
and 04. 
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FIG. 5: The T dependence of the total observed torque mag- 
netization Aloba measured at _ff = 5 T in sample 06 showing 
the weak T-linear change above Tonset (arrows), the onset 
temperature of diamagnetism. The linear extrapolation (dot- 
ted lines) below Tonaet is used to determine the anisotropic 
van Vleck contribution Ax"^*"!!. Below Tonset, the magnitude 
\Md\ increases rapidly. 




FIG. 6: The curves of M'^^i, (observed magnetization with van Vleck term removed) in samples 03, 04, 05, 055 and 06 (Panels 
a, b, c, d and e, respectively). Each of these curves is comprised of a diamagnetic term Md{T,H) that is dominant in low 
H (<10 T), but is suppressed above 20-30 T, and a paramagnetic spin term l\Ms{T, H) that is nearly constant in H in high 
fields. As we increase x, the diamagnetic term Md becomes increasingly dominant and harder to suppress in high fields (note 
expanded field scale in Panel d). Samples 055 and 06 show broad Meissner transitions with Tc ~ 0.5 and 5 K, respectively. 
The other samples show no Meissner transitions. For clarity, we have not displayed curves at T below 4 K (see Fig. 2 in [|). 
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FIG. 7: Curves of the the diamagnetic term Md{T,H) and 
the paramagnetic spin term AMs (T, H) in sample 03 {x = 
0.03) at fixed T as indicated. At each T, after removing 
the orbital van-Vleck term, the remaining signal is given by 
M'^^^{T,H) ^ AM,{T,H) + Ma{T,H). The plotted curves of 
Md and AMs add to reproduce M'^f^^. The inset shows the 
T dependence of the prefactor ViT) obtained from the fit to 
Eq. 13 
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FIG. 8: Comparison of the hysteretic loops in Samples 06 (Panel a) and 055 (b). In spite of the small difference in x, the 
hysteresis loops are much larger in 06 than in 055 measured at the same T. This reflects the rapid decrease of the zero-Kelvin 
irreversibility field Hirr{0) as x approaches Xc ~ 0.055. The inset (Panel c) shows an expanded view of the hysteresis at 0.35, 
0.7 and 1.5 K in Sample 055. 
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FIG. 9: Jumps in the field sweep-up portion of the hysteretic 
curves of Md vs. H (arrow at Havai) in sample 06 at T = 
0.35 K at a field sweep rate of 5 T/min. The jumps are not 
observed when the sweep rates are reduced by a factor of 10. 
Vortex avalanches are observed only in the vortex-solid phase. 



